Non-Oberbeck-Boussinesq effects in turbulent thermal convection in ethane close to 

the critical point 
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As shown in earlier work (Ahlers et at, J. Fluid Mech. 569, 409 (2006)), non-Oberbeck Boussinesq 
(NOB) corrections to the center temperature in turbulent Rayleigh-Benard convection in water and 
also in glycerol are governed by the temperature dependences of the kinematic viscosity and the 
thermal diffusion coefficient. If the working fluid is ethane close to the critical point the origin 
of non-Oberbeck-Boussinesq corrections is very different, as will be shown in the present paper. 
Namely, the main origin of NOB corrections then lies in the strong temperature dependence of the 
isobaric thermal expansion coefficient /3(T). More precisely, it is the nonlinear T-dependence of the 
density p(T) in the buoyancy force which causes another type of NOB effect. We demonstrate that 
through a combination of experimental, numerical, and theoretical work, the latter in the framework 
of the extended Prandtl-Blasius boundary layer theory developed by Ahlers et at, J. Fluid Mech. 
569, 409 (2006). The latter comes to its limits, if the temperature dependence of the thermal 
expension coefficient /3(T) is significant. The new measurements reported here cover the ranges 
2.f < Pr < 3.9 and 5 x f0 9 < Ra <2x f0 12 and are for cylindrical samples of aspect ratios f .0 and 
0.5. 



I. INTRODUCTION 

Fluid motion in the presence of temperature gradi- 
ents is an important phenomenon in nature and indus- 
trial processes. Among the many examples are oceanic 
streams, cloud motions, and gusts of wind that can be 
felt on a human scale. The paradigmatical system for 
such thermally driven flows is the Rayleigh-Benard (RB) 
setup: a fluid- filled container heated from below and 
cooled from above, fn this classical problem, the flow 
is determined by the scale and geometry of the con- 
tainer, the material properties of the working fluid, and 
the top-down temperature difference A = T\, — T t > 0. 
In the last two decades, considerable progress has been 
achieved in our understanding of global and local prop- 
erties and the flow organization of turbulent RB con- 
vection, through a combination of experimental (see e.g . 

h, a a a ib a a m m el m, m m ee m 

numerical (see e.g. 
H, HI), and theo- 
retical work (see e.g 

The temperature difference in a RB cell can be in- 
creased in a controlled way. However, in principle the 
transport coefficients of the fluid can depend on the local 
temperature and density and thus vary across the height 
L of the container. Since space-dependent properties of 
such kind are undesirable in first instance, one tends to 
restrict the convection regime to sufficiently small inter- 
vals of A. But even so, further simplifications are pro- 
gressively required in the analysis of RB convection. In 
this spirit, a standard approximation due to Oberbeck 
[39l | and Boussinesq assumes that (see also [4ll. l42j ) : 



OB.l The dynamic viscosity 77, the thermal conductiv- 
ity A, the thermal expansivity /?, and the isobaric 
specific heat cp are constant throughout the fluid. 

OB. 2 Density variations are taken into account only in 
the buoyancy force term. 

OB. 3 The temperature dependence of the density p is 
linearized in the buoyancy force as: 



P(T) 
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where T m = (X& + T t )/2 is the arithmetic mean 
temperature between the plates and X m = X{T m ) 
denotes the fluid property X evaluated at T m . 

Next to the aspect ratio, within the OB approximation 
two dimensionless parameters characterize the RB flow: 
The Prandtl number Pr = v m j ' K m follows from the ra- 
tio between the kinematic viscosity v m = r\ m j p m and the 
thermal diffusivity K m = A m /(/9 m cp iTn ,). The dimension- 
less thermal driving can be conveniently represented by 
the Rayleigh number Ra = /3 m gL 3 A/ (i/ m K m ) , where g 
denotes the gravitational acceleration. 

The manner in which high Rayleigh numbers are 
achieved is crucial for the emergence of non-Oberbeck- 
Boussinesq effects (NOB). Since turbulent convection 
may involve spatiotemporal changes in the fluid proper- 
ties, considerable efforts have been devoted to the identi- 
fication of dominating sources of NOB effects. In liquids 
like water [43| and glycerol [4^ ]. (45|, for example, NOB 
effects are dominated by deviations from (OB.l) since 
the viscosity strongly decreases with temperature. On 
the other hand, when the working fluid is gaseous ethane 
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[49| . deviations from (OB.l) and (OB. 2) lead to NOB 
effects stronger than those in the aforementioned liquids. 

In the present study, we shall focus on deviations from 
(OB. 3) by considering the nonlinear temperature depen- 
dence of the buoyancy force. In particular, ethane close 
to its critical point [47 ] is chosen as the working fluid 
and the temperature T c in the center of the container is 
measured as indicator of NOB effects. 

There are two possibilities to characterize the physics 
beyond condition OB. 3 as described by eq. {1}, if the den- 
sity p(T) has a strong nonlinear T-dependence. First, one 
can introduce a T-dependent thermal expansion function 
P(T) instead of /3 m , defined in terms of the density p(T) 
by 



p{T) = p m - Pm(T - T m )/3(T). 



(2) 



Second, one refers to the common isobaric thermal ex- 
pansion coefficient /3(T), defined as usual by 



f3(T) 



1 dp(T) 
p(T) dT 1 



(3) 



which now is temperature dependent. Both are related 
by (3{T) = -d log[l — (T — T m )(3{T)\ / dT. Under the 
condition OB. 3 of linear T-dcpcndcncc of p the thermal 
expansion function is constant, (3 = /3 m , while the ex- 
pansion coefficient f3(T) is given by (3(T) — /3 m /[l — (T — 
T m )/3 m ], still depending on temperature. Of course all 
three coincide at T m . An advantage of considering the 
thermal expansion coefficient f3(T) is that it is a well de- 
fined thermodynamic derivative. The advantage of /3(T) 
on the other hand is that it immediately reflects the non- 
linear T-dependence of p(T). Also f3(T) refers to a single 
thermodynamic state and describes the local T-slope on 
an isobar (normalized by the local density), while /3(T) 
refers to a pair of states, namely to the reference state 
T m in addition to T and describes the secant to the p{T)- 
curve (normalized by the reference density). /3(T) will 
therefore in general vary less with T than f3(T). 

It will turn out that it is the significantly different In- 
dependence of $(T) (or f3(T)) on the two sides of the 
critical isochore of ethane, which leads to opposite shifts 
of the center (bulk) temperature T c , yielding T c < T m on 
the gas-like (i.e. high-temperature, see Fig.Q]below) side 
and T c > T m on the liquid-like (low-temperature) side. 
On the gas-like side $(T) increases from bottom to top 
and on the liquid-like side it decreases. 

Our approach consists of three stages: boundary layer 
(BL) theory, experiments, and direct numerical simula- 
tions (DNS). First, we address in section [TT] an extension 
of boundary layer theory that considers deviations from 
(OB.l) and (OB. 2). Even though the buoyancy force 
is not included in the BL equations (only the longitu- 
dinal momentum is taken into account here), we com- 
pute T C (A) for several pressures P m . Then, experimen- 
tal measurements of T C (A) are presented in section UlI Bl 
and compared with BL results in section HVl Given the 



significant discrepancies between part of them, we ad- 
dress in section [V] direct numerical simulations that ex- 
plicitly consider deviations from (OB. 3). In particular, 
for T m = 27° C and P m = 51.72 bar, it is shown that 
NOB effects in ethane are dominated by the nonlinear 
dependence of the buoyancy force on temperature. Fi- 
nally, our conclusions are summarized in section PVTl Ap- 
pendices [5] and [B] are devoted to the derivation of the 
boundary layer equations with variable transport coeffi- 
cients and Appendix [U] compiles the Nusselt number cor- 
rections for the real and various hypothetical ethane-like 
fluids in a table. 



II. BOUNDARY-LAYER THEORY 

A central aspect in Rayleigh-Benard convection is the 
understanding of the boundary layers formed along the 
bottom and top plates. Though they preserve a lami- 
nar character for Ra < 10 12 , their instabilities impact 
the Nusselt number Nu (the effective heat flux relative 
to thermal conduction A m A/L) |H,|4^|. As reported in 
references [H, 0, Hgj], BL flows of this nature are sig- 
nificantly influenced by the coupling between the fluid 
properties and the temperature gradient across the con- 
tainer. In particular, it was shown that NOB effects on 
T c can be reasonably described by extending the Prandtl- 
Blasius boundary layer theory [Hoi . I5l1 ] . Next we review 
such an extension and further assess its intrinsic limita- 
tions. 

Assume that the density p, the temperature T, and 
the velocity u are stationary fields, which depend only on 
the longitudinal x and transverse z coordinates. Then, 
under the boundary-layer approximation, we write the 
continuity and the ^-momentum equations as (see also 
appendix IA ip : 



^{pu x } + ±{pa M } 



= 0, 



du x 
dx 



du x 
dz 







d [ du x 
V 



dz 



(4) 
(5) 



Here z measures the vertical distance from the bottom or 
top plates, respectively, the velocity components at z — 
are subject to no-slip boundary conditions: u x (x, 0) = 
and u z (x, 0) = 0. Moreover, in the bulk of the flow, u x is 
supposed to match the large scale wind velocity U c in the 
center (bulk) of the RB sample [H^], i.e., u x (x, oo) = U c . 
Note that within the BL theory we cannot calculate U c ; 
here we only have to assume that it is the same close to 
the top and the bottom BL which is supported by our 
numerical simulations reported in section V. 

In the same spirit, the temperature field T(x, z) is gov- 
erned by (cf. appendix IA 2p : 
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dT 

dz 



du x du z \ 

a ]^ + ^' (6) 
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where 7 = cpjcy is the ratio between the isobaric 
and isochoric specific heats and (3 = —p~ 1 dp/dT de- 
notes the isobaric thermal expansion coefficient. At the 
plates T(x,0) = T&t and in the bulk (center) of the flow 
T(x, 00) = T c . 

The coupling between the bottom and top boundary 
layers is determined (cf. (43[) by the heat fluxes Qb,t 
through the plates, considered to be equal: 



Qb 



-A b 



dT 

dz 



= -At 



dT 

dz 



= Qt 



This condition establishes an implicit dependence of the 
center temperature T c on the heat fluxes Qb = Qt = Q- 
Note again that both the dynamic viscosity 77 and the 
heat conductivity A depend on both temperature and 
density i.e., r](T,p) and A(T, p). Before addressing the 
technicalities around the (numerical) integration of equa- 
tions ([I])-©, we shall benefit from a key argument in 
boundary-layer theory: Prandtl's self-similar ansatz. 



Here fj = rj/rjm is the dimensionless viscosity, whose in- 
dependence 77 ' is given by 



V' 



*n e'+flf] p' 



dO 



dp 



where = (T — T t )j A denotes the dimensionless tem- 
perature. 

Next, assuming that the pressure P m is constant 
throughout the fluid, one finds 



-p/30', 



(12) 



with (3 = (3 A. The boundary conditions at the respective 
walls are p(0) = pb,t, p'(0) = - p b ,t Pb,t ' b ,v and 
p(oo) = p c . 

Finally, we also write the temperature equation ([6]) in 
self-similar form as (see Appendix IB 2p 

A6" + j^cp Prf + A'j 6' = 0, (13) 



A. Self-similarity 

Because of the y-independence, assumed in the Prandtl 
BL theory the boundary layer flow is mathematically a 
2D flow. Therefore the system of partial differential equa- 
tions <(4j)— ((6j) for the BL flow can be reduced to ordinary 
differential equations (ODEs) by introducing a stream 
function We do this differently from the usual proce- 
dure by including in its definition the density in order to 
automatically fulfil the continuity equation by construc- 
tion. 



pu x 



pu z = 



a* 

dz ' 
dx ' 



(8) 
(9) 



where p = p/p m is the density nondimensionalized with 
Pm = p{T m , P m ). Apparently the continuity equation au- 
tomatically follows from ©^(O- Next, cf. Appendix [Bl 
we may introduce a self-similarity variable Z = z/i c (x) 
and a similarity function ^>(Z) = ^(x, z)/(£ c (x) U c ), 
such that £ c (x) = y/ xv m /U c . Thus the velocity com- 
ponents are 



P 




(10) 



with boundary conditions ^(0) = = ^'(O) and 
*'(oo) =p c . 

In terms of (fit))) , the viscous BL equation (O can be 
written as: 







2 p 



1L* 
2 p 



p" 



^'U(ii) 



where A = A/A m and cp = cp/cp, m . Equation (|13p is 
subject to 6(0) = &b,t and 9 (00) = C . 



B. Results 

The coupled ODEs (n])-CL3|) with the respective 
boundary conditions and the heat-flux conservation ([7|) 
are solved numerically with a shooting method [53[ . The 
integration domain is restricted to A-intervals where the 
transport properties are concave/convex functions of the 
temperature. In particular, we have chosen ethane as 
the working fluid since its properties are known very well 
[4^ . even close to its critical point (T*, P», p*) [see figure 
[T]. All material properties 77, A,p, /3, and c p are imple- 
mented in their full dependence on T In this manner, the 
computation of temperature and density profiles does not 
involve any fit parameter. 



1. Vertical profiles 

An insight into the structure of the BLs can be 
achieved by studying typical profiles along the z- 
direction. To describe them, let us consider a representa- 
tive case in which the pressure is fixed at P m = 0.849 -P*, 
the mean temperature at T m = 40 °C, and the thermal 
difference between the plates at A = 15 K. 

In figure [21 the temperature and density p are plot- 
ted as functions of the similarity variable Z . As shown 
in panel |2^,, the center temperature C is smaller than 
the mean temperature m = 0.5, clearly indicating a 
top-down symmetry breaking. Such symmetry breaking 
is also reflected in the density profiles shown in panel 
since the center density p c is larger than the mean 
density p m — 1. Notwithstanding the pronounced curva- 
tures in Q(Z) and p(Z), we shall restrict our attention 



4 



(a) 





T/T 



1 2 3 4 5 6 



FIG. 1: Pressure-temperature plane in reduced units. Star: 
critical point of ethane (T„ = 32.18°C, P* = 48.718 bar). 
Heavy line: liquid-vapor coexistence curve. Dotted line: crit- 
ical isochore. The horizontal arrows show the maximum 
temperature intervals A, centered at T m = 27 °C (left) and 
T m = 40 °C (right). The pressures are P m /P, = 0.849, 0.991, 
1.026, 1.062, and 1.104 (bottom to top). 



(b) 



to the asymptotic value O c 
NOB effects. 



as a convenient indicator of 



2. Center temperature 

To compute the difference T c — T m as a function of A, 
we have chosen a particular set of isobars in the phase 
diagram of ethane. As shown in figure [TJ our selection 
of A-intervals falls into two classes: (i) those intervals 
centered at T m = 27 °C and (ii) those centered at T m = 
40 °C. 

As for the latter (the more gaseous case) , figure 
shows that the center temperature is a decreasing func- 
tion of A. The top-down symmetry of the flow is bro- 
ken in such a way that the top boundary layer tends to 
become thinner than its bottom counterpart, eventually 
leading to a temperature reduction in the center of the 
flow. Though this result has been originally reported and 
explained already in reference (46j , we briefly mention it 
here for completeness of discussion. 

Focussing now on the class of A-intervals centered at 
T m = 27° C, figure|4jD shows that the center temperature 
becomes larger than the mean temperature between the 
plates. Such NOB effect is different to what we found in 
reference (46j and repeated in figure 2] a, where we have 
focused on the more gas-like case. 

To understand this we argue that the intervals under 
consideration (centered at T m — 27° C) now correspond 
to a region of the phase diagram, where the material 
properties of ethane behave more similar to those of the 
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Bottom BL 
Top BL 



ICL 




FIG. 2: (a) Temperature and (b) density profiles at P m = 
0.849 ■ P., T m = 40° C and A = 15 K. The t hermal sl ope 
thicknesses at the bottom or top \f it /L — ab,t/\/ v m /LU C are 
seen to have prefactors of about « 2.8 and at ~ 2.5, i.e., 

A si \ \Sl 
b > A i • 



liquid phase. NOB effects in classical liquids (such as 
water and glycerol) were already discussed in references 
[43L H, 45] . One of our aims in the present work is to fur- 
ther assess the differences in the NOB effects between the 
more liquid-like versus the more gas-like fluids, see Sec. 
IIII1 It will turn out that this will show us the limitations 
of boundary-layer theory, see Sec. HVl To this end, we 
first consider now additional experimental details on the 
A dependence of the center temperature T C (A). 



III. EXPERIMENT 
A. Apparatus 

The apparatus was described in detail before in 
Ref. [54j |. where a schematic diagram is shown in Fig. 
2 . Here we give a brief description and details specific 
to the present high-pressure sample cell shown in Fig. [3] 
Working from the inside out, the sample cell was sur- 
rounded by a can containing ambient air. The air space 
inside the can was filled with low-density open-pore foam 
to prevent convection outside the sample. The maximum 
possible diameter of the sample top plate was 10 cm, al- 
lowing for inside sample diameters typically up to about 
D = 8 cm. The entire apparatus was of sufficient length 
to accommodate a sample with L ~ 16 cm (r ~ 0.5). 
Heat was applied at the sample bottom by a metal-film 
heater covering the entire active bottom-plate area uni- 
formly. 

The top plate was cooled by a circulating water bath. 
The water was cooled when passing through a heat ex- 
changer external to the main apparatus which in turn 
was cooled by a separate water circuit driven by a Nes- 
lab or Lauda refrigerated circulator with a temperature 
stability of 0.01°C. Just before entering the apparatus, 
the water was heated by a heater consisting of about 25 
m of teflon-insulated AWG30 (0.5 mm dia) copper wire, 
stuffed into the inlet line and thus immersed in the water. 
The large contact area between the water and the heater 
wire provided excellent heat exchange and uniform heat- 
ing of the water. The heater was computer controlled in 
a feedback loop with a thermometer located in the top 
plate of the sample cell. The bath-temperature stability 
achieved in this way was a few tenths of a milli-Kelvin. 
The water entered the bottom center of the apparatus, 
flowed upward through an annular channel around the 
can, and was distributed over the top plate by a set of 
jets. With this arrangement the entire can was kept at 
the top-plate temperature and parasitic heat loss from 
the side wall and the bottom plate due to conduction 
through the air/foam as well as by radiation was reduced 
to a level that was negligible compared to the heat trans- 
port by the convecting fluid. 

After cooling the top plate the water returned through 
an annular channel located just outside of and mildly 
insulated from the incoming channel. Since the water, 
while cooling the top plate, was never heated by more 
than a few mK, the returning water provided an excellent 
adiabatic thermal shield at the top-plate temperature, 
thus stabilizing the interior temperatures and preventing 
significant variations in time of the parasitic heat losses 
from the bottom plate. The entire apparatus sat on a 
chlorinated poly (vinyl chloride) (CPVC) base plate with 
appropriate channels and feed-throughs to accommodate 
the water circuit, the electrical leads, and the fill capillary 
going to the sample. 

One of two high-pressure sample-cells (see Fig. [3]) was 
installed inside the can. It had the shape of a cylin- 
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FIG. 3: Schematic diagram (approximately to scale for F — 
1) of the high-pressure sample cell surrounded by the can 
containing ambient air and foam. 



der with D = 7.63 cm. One cell had an internal 
length L = 7.62 cm, corresponding to an aspect ratio 
r = D/L = 1.00. Another one had L = 15.24 cm, yield- 
ing r = 0.500. The top and bottom plates consisted of 
thick copper disks. Each of the two plates had an anvil, 
1.59 cm thick, of diameter essentially equal to D. that 
was a close slide fit in the type 4340 steel side wall. After 
machining, the side wall was heat treated at 830°C and 
oil quenched. This procedure is expected to lead to a ten- 
sile strength of about 13 kbars. The side wall had a flange 
at each end, of thickness 0.95 cm. Each flange was bolted 
(not shown in the figure) and "O" -ring-sealed to one of 
the copper end plates. A top and bottom thin section of 
the side wall had a thickness of 0.051 (0.076) cm for the 
r = 1 (r = 0.5) cell. This thin section overlapped the 
copper anvils and extended into the sample region by 0.95 
cm. Connecting the thin section was a central section of 
wall thickness 0.15 cm that provided enhanced strength; 
since the turbulent system contained only a very small 
thermal gradient in its center, the thicker wall section did 
not significantly enhance the wall heat-transport. The 
sample entered the bottom copper plate through a cap- 
illary from the side, and then proceeded through a very 
small hole (shown in white on the left side of the figure) 
into the gap between the bottom-plate anvil and the side 
wall. 

The sample was connected to a manifold through a 
capillary. Also connected to the manifold was a separate 
pressure-regulation volume of 600 (1000) cm 3 in the case 
of the r = 1 (r = 0.5) cell that could be heated above 
the ambient temperature by a heater wrapped around its 
outside. The temperature of this "hot volume" was con- 
trolled in a feedback loop with a pressure gage [55[ . The 
pressure stability typically was better that one milli-bar. 
The entire system was designed to safely withstand pres- 
sures up to 60 bars. For the pressure measurements we 
used a Paroscientific model 745 pressure standard with 
an accuracy or 80 ppm (about 6 milli-bars) and a reso- 
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lution of 1 ppm (about 70 micro-bars) of full scale. 

A substantial fraction of the heat current passed 
through the side wall. This current was measured for the 
evacuated cell and subtracted from all other measure- 
ments; but as was recognized some time ago ([IB, l57j). 
this is not an adequate procedure because of the height- 
dependent temperature gradients that prevail in the wall 
when the cell contains turbulently convecting fluid. We 
did not attempt a correction for this non-linear side-wall 
effect in the present case because we do not believe that a 
reliable correction is possible when the side-wall conduc- 
tance is large. For this reason our values of Nu under OB 
condition are about 25% larger than other measurements 
at similar Prandtl numbers ([16]). However, we believe 
that the deviations of Nu and of T c from their Boussi- 
nesq values, which (as we shall see below in Sect. IIII B 2|) 
depend primarily on the nature of the top and bottom 
boundary layers rather than on the fluid interior, were 
obtained reliably. 

The top and bottom temperatures T t and T b were 
determined from the average of six thermistors imbed- 
ded close to the fluid in each of the top and bottom 
plates ([El]). These thermistors were calibrated against 
a platinum-resistance thermometer purchased from Hart 
Scientific. This thermometer was supplied with a cali- 
bration, accurate to ±7 mK, on the ITS90 temperature 
scale. The average temperature readings were used to 
obtain A = T b - T t and T m = (T t + T b )/2. Small correc- 
tions for the temperature gradients in the copper plates 
were applied. The center temperature T c was taken to 
be the average of the temperatures measured with eight 
thermistors attached to the outside of the side wall at the 
horizontal mid-plane, equally spaced in the azimuthal di- 
rection (see, for instance, Ref. [HJ or [Bl^]). 

In order to obtain an estimate of the OB values of the 
Nusselt numbers, a power law Nuob — N^R^'fi was 
fitted to the Nusselt-number measurements at relatively 
small A where /3 m A < 0.05, adjusting ~f e ff and N . Such 
fits yielded values of 7 e // close to 0.30. All the measured 
values of Nu regardless of (3 m A were then divided by the 
power-law value at the measured Raylcigh numbers to 
give Nu/ Nuob at all A. 

All measurements reported here were made with many 
values of A at each of a few constant values of T m and 
P. In both cells we used ethane at elevated pressures as 
the fluid. The thcrmophysical properties were calculated 
from the formulas given in Ref. |47| . For extensive dis- 
cussions of the uncertainties of these properties we refer 
to that paper. It is difficult to determine the absolute 
errors for the Rayleigh and Nusselt numbers that results 
from property uncertainties, but we expect that an es- 
timate of a few percent is not unreasonable. Since in 
the present paper we are concerned only with the ratios 
Nu/ Nuob, and since all data are taken as a function of 
A at a given mean temperature and pressure and evalu- 
ated at the same P and T m , property errors cancel to a 
very large extent. 
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FIG. 4: Deviation T c — T m as function of A, for (a) T„, 
40 °C and (b) T m = 27 °C. 



B. Results 

1. The Nusselt number Nu(Ra) 

As indicated above, we do not regard the results for 
Nu(Ra) to be very accurate because of unknown effects 
due to the relatively large wall conductivity. Nonetheless 
we show the results for T = 1 at several T m and P in 
Fig. on logarithmic scales. Over a wide range of Ra 
one sees that they are a few percent higher than the re- 
sults from Refs. [6[ and 0], and we attribute this to the 
influence of the side-wall conductivity on our data. At 
the largest Ra our results increase more rapidly with Ra, 
and data at different T m and P begin to differ from each 
other. We attribute this phenomenon to NOB effects. 
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FIG. 5: (a): The Nusselt number Nu as a function of the 
Rayleigh number Ra for T — 1.00. Solid circles: P = 51.72 
bars and T m = 40°C (Pr = 2.58). Open circles: P = 51.72 
bars and T m = 27°C (Pr = 2.99). Solid squares: P = 51.72 
bars and T m = 24°C (Pr = 2.71). Open squares: P = 51.72 
bars and T m = 31° C (Pr = 3.85). Solid diamonds: P = 53.79 
bars and T m = 40° C (Pr = 3.79). Open diamonds: P = 50.00 
bars and T m = 40° C (Pr = 2.09). Stars: from Ref. for 
r = 0.5 after correction for sidewall effects. Plusses: from 
Ref. for P = 0.5. 
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FIG. 7: The ratios X/X m for P = 51.72 bar (P/P, = 
1.062) of several property values X at temperatures T — T m 
to the value of X at T m (based on Ref. [47j). (a) and 
(c): T m = 27.00° C. (b) and (d): T m = 40.00° C. (a) and 

(b) : thermal conductivity A (short dashed line), density p 
(dotted line), and dynamic viscosity r\ (dash-dotted line). 

(c) and (d): thermal expansion coefficient f3 (solid line), 
j3 (long dashed line), and heat capacity at constant pres- 
sure cp (double-dashed dotted line). The reference values 
X m for 27(40)°C are A m = 0.07328(0.04343) W m _1 K _ \ 
p m = 331.12(123.26) kg m" 3 , rj m = 4.030(1.502)10'° kg s" 1 
m" 1 , p m = 0.01649(0.03815) K" 1 , c P , m = 5434(7452) J kg" 1 
K -1 and mean Prandtl number Pr m — 2.99(2.58). 
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2. Aspect-ratio dependence 

In Fig. [6] we compare results obtained at a mean tem- 
perature T m = 40.00°C and pressure P = 51.72 bar 
(P/P* = 1.062, Prandtl number Pr = 2.58) in the sam- 
ple of aspect ratio F = 0.50 (open circles) with previously 
reported results ([11]) for T = 1.00 (solid circles). One 
sees that the NOB effect on Nu and on the center tem- 
perature T c is within our resolution independent of T. 
This shows, as expected, that the NOB effects are con- 
fined essentially to the boundary layers. The length of 
the sample interior, which is nearly isothermal (see, how- 
ever, Ref. [59]) regardless of its length, does not have a 
large influence. 



FIG. 6: (a): The ratio of the measured Nusselt number Nu to 
the estimate Nuob of the Nusselt number under Boussinesq 
conditions as a function of the applied temperature difference 
A. (b) The deviation of the center temperature T c from the 
mean temperature T m as a function of A. All measurements 
were made at T m = 40.00° C and a pressure of 51.72 bar 
{P/P, = 1.062) where the Prandtl number is 2.58. Open 
symbols: V = 0.50. Solid symbols: data from Ref. [4(| with 

r = i.oo. 



3. Dependence on fluid properties 

Interesting insight into the influence of various prop- 
erty variations with temperature can be gained by mea- 
suring T c and Nu along an isobar on the two sides of 
the temperature T^(P) at which the critical isochore is 
reached. In Fig. [7] we show the variation along the isobar 
P = 51.72 bar = 1.062 P* of the thermal conductivity A, 



8 




0.4 



A (K) 



FIG. 8: (a): The ratio of the measured Nusselt number Nu 
to the estimate Nuob of the Nusselt number under Boussi- 
nesq conditions as a function of the applied temperature dif- 
ference A. (b) The deviation of the center temperature T c 
from the mean temperature T m . These measurements were 
made for F = 1.00 at a pressure of 51.72 bar (P/P* = 1.062). 
Open symbols: T m = 27.00° C (Pr = 2.99). Solid symbols: 
T m = 40.00° C (Pr = 2.58). These two temperatures are on 
opposite sides of, but not equidistant from, the temperature 
T$ = 34.97° C where the critical isochore is reached on this 
isobar. 



density p, dynamic viscosity r/, thermal expansion coeffi- 
cient (3, and heat capacity cp for the cases T m — 27.00°C 
(left panels) and 40.00°C (right panels). These two val- 
ues are on opposite sides of but not quite equi-distant 
from = 34.97°C. In the upper two panels one sees that 
the variations of A, p, and r\ are relatively small, have the 
same trends with T though quantitatively they are some- 
what different on the two sides, with maximum changes 
by less than a factor of two over temperature ranges that 
are small enough to avoid including T^. On the other 
hand, the expansion coefficient and the heat capacity 
(lower two panels) vary by a factor of five or more. Thus 
one might expect them to dominate the NOB effects. In- 
terestingly they have opposite trends with T — T m , the 
temperature derivatives of both (3 and cp are positive 
below and negative above T^; that is below the criti- 
cal isochore (along the temperature axis), on the more 
liquid-like side, (3 and cp are smaller at the top (colder) 
than at the bottom (warmer) end of the sample, with this 
relationship reversed above the critical isochore, on the 
more gas-like side, where (3 and cp decrease from bottom 
to top. 

In Fig. [5] experimental results are presented for T = 1.00 
at a pressure P = 51.72 bar = 1.062 P* . They are for 
the two mean temperatures T m — 27.00 (open circles) 
and 40.00°C (solid circles) of Fig. where the Prandtl 
numbers are 2.99 and 2.58 respectively. In both cases 
we used A values sufficiently small so that T t (T&) only 
reaches down (up) to so that the applied temperature 
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FIG. 9: The ratios X/X m for P = 55.17 bar (P/P, = 
1.132) of several property values X at temperatures T — T m 
to the value of X at T m (based on Ref. [47j). (a) and 
(c): T m = 35.00°C. (b) and (d): T m = 41.00°C. (a) 
and (b): thermal conductivity A (dashed line), density p 
(dotted line), and dynamic viscosity r\ (dash-dotted line), 
(c) and (d): thermal expansion coefficient /3 (solid line), 
j3 (long dashed line), and heat capacity at constant pres- 
sure cp (double-dashed dotted line). The reference values 
X m for 35(41)°C are A m = 0.06674(0.05163) W m" 1 K~\ 
p m = 282.48(153.43) kg m" 3 r) m = 3.168(1. 726)10~ 5 kg s^ 1 
m" 1 , m = 0.04177(0.06876) K _1 , c P<m = 9617(12534) J 
kg" 1 K" 1 , and mean Prandtl number Pr m — 4.56(4.20). 



difference does not straddle T$. One sees that the NOB 
effects increase Nu on both sides of the critical isochore. 
On the high-temperature side (solid circles) the NOB ef- 
fect is larger for the same A. This is consistent with the 
larger variation of the fluid properties at equal values of 
T — T m revealed above in Fig. [7] 

The NOB effect on T c is of opposite sign on the two sides 
of the critical isochore. For T < (open circles) NOB 
conditions increase T c above T m , whereas for T > 
(solid circles) T c is reduced below T m . This observation, 
in conjunction with the properties shown in Fig. sug- 
gests that for these fluids the temperature drops A tj t, 
across the boundary layers are determined primarily by 
(3 and/or cp, with A t > Ah (A t < A&) when (3 and/or 
cp are smaller (larger) at the cold top end of the sample 
then at the warm bottom end. As pointed out before, [431 ] 
for the Nusselt number the contributions to the thermal 
resistance at the two boundary layers add, and it does 
not matter much whether the larger or smaller contribu- 
tion comes from one end or the other. Thus, for Nu the 
NOB effect is in the same direction in both cases. As 
was the case for Nu, the NOB effect revealed by T c is 
larger in magnitude above than it is below. Again we 
attribute this difference primarily to the difference in the 
variations of the properties shown in Fig. [7] 
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FIG. 10: (a): The ratio of the measured Nusselt number 
Nu to the estimate Nuob of the Nusselt number under 
Boussinesq conditions as a function of the applied temper- 
ature difference A. (b) The deviation of the center tem- 
perature T c from the mean temperature T m . These mea- 
surements were made for T = 0.50 at a pressure of 55.17 
bar (P/P* = 1.132). Open symbols: The mean tempera- 
ture T m = 35° C (Pr = 4.56).' Solid symbols: T m = 41° C 
(Pr = 4.20). These two temperatures are on opposite sides 
of and nearly equidistant from the temperature T^, = 38.06° C 
where the critical isochore is reached on this isobar. 



IV. COMPARISON BETWEEN 
BOUNDARY-LAYER THEORY AND THE 
EXPERIMENTAL RESULTS 

We now compare the experimental measurements from 
the previous section with the boundary layer results pre- 
sented in section [iTJ In particular, since the compar- 
ison with experiments at T m — 40° C (more gas-like 
ethane) was already discussed in reference [46j], finding 
good agreement between experiment and the extended 
BL theory, we devote special attention to the new mea- 
surements at T m = 27° C and P/P* = 1.062 (more liquid- 
like ethane). 

As shown in figure 111) the curve for more liquid-like 
ethane obtained from BL theory considerably deviates 
from the experimental data. This is remarkably different 
from the comparison of BL theory with gaseous ethane, 
presented in [461 ] . where instead a good agreement was 
observed. This suggests that even though BL theory rea- 
sonably captures NOB effects associated with (OB.l) and 
(OB. 2), further corrections are essential in the present 
liquid-like ethane case. Among them, deviations from 
(OB. 3) seem to be the natural candidate for the failure 
of BL theory, since the buoyancy force is not included in 
the BL equations but apparently affects the thermal con- 
vection. Thus, in order to reveal the importance of non- 
linear buoyancy in thermal convection, we shall perform 
direct numerical simulations (DNS) of the RB problem. 



In Fig. [9] we show the variation along the isobar P = 
55.17 bar (P/P* = 1.132) of the various properties for 
the cases T m = 35.00°C (left panels) and 41.00°C (right 
panels). These two temperatures are also on opposite 
sides of and nearly equi-distant from the critical isochore, 
for this pressure at = 38.06°C. Again the variation of 
the expansion coefficient and the heat capacity (lower two 
panels) is much larger than that of the other properties. 
At a given \T ~ T m \, all the variations are more similar 
in magnitude on the two sides of than they were for 
the case of Fig. [7] 

In Fig. [10] experimental results corresponding to the 
conditions of Fig.[9]are presented for T = 0.50. They are 
for the two mean temperatures T m — 35.00 (open circles) 
and 41.00°C (solid circles) where the Prandtl numbers 
are 4.56 and 4.20 respectively. In both cases we used 
A < 6 K so that T t (Tb) reaches down (up) to when 
T m = 41°C (35°C) while the applied temperature differ- 
ence does not straddle T^. For this case one sees that 
the NOB effects on Nu are similar on the two sides of 
the critical isochore. Again, the NOB effect on T c , al- 
though of about the same magnitude, is of opposite sign 
on the two sides. For T < (open circles, more liquid- 
like) NOB conditions increase T c above T m , whereas for 
T > T<f, (solid circles, more gas-like) T c is reduced below 



V. DIRECT NUMERICAL SIMULATIONS 

As shown in refs. [45l [60]|. two-dimensional direct nu- 
merical simulations may be useful for the study of the 
tiny NOB effects which occur in RB convection in liq- 
uids. In particular, even with the restrictions to two- 
dimensional geometry and to incompressibility of the 
fluid flow, the effects on the center temperature and the 
Nusselt number shift could be reasonably captured in 
the cases of water and glycerol. However, the liquid- 
like ethane just above the critical pressure has a stronger 
temperature dependence of the density than water and 
glycerol. To quantify this, a comparison between the 
ethane properties around T m = 27°C and P/P* = 1.062 
with water and glycerol around T m = 40° C is reported in 
figure fL2l For the case of ethane the incompressible flow 
approximation seems to be questionable, or at least less 
justified. But we will show that a dopting the same ap- 
proach used for water and glycerol [43, [6(| also proves to 
be useful to study NOB effects in ethane and the results 
are well consistent with experiment. 

Further insight into the problem is given by considering 
several cases of artificial ethane-like fluids, namely fluids 
which have only one, or some, of their material prop- 
erties dependent on temperature, while the others are 
kept constant. In particular, as discussed in the previous 
section, we will examine the relevance of the nonlinear 
temperature dependence of buoyancy on the center tern- 
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FIG. 11: Deviation T c - T m as function of A, for T m = 27 °C 
and P/P, = 1.062. The symbols (x) correspond to exper- 
imental measurements. The dotted line is obtained from 
boundary layer theory. The red symbols (O, A, and □) 
with error bars correspond to the incompressible DNS results 
described in section [V] measured at different Ra — 10 6 to 
Ra = 10 s . Note that though these Rayleigh numbers are 
smaller than in the experiments (Ra = O(10 9 )-O(10 10 )), the 
comparison is still appropriate because the T c shift has proven 
to be rather independent of Ra for given A, provided one is 
beyond the onset of the chaotic motion at Ra ~ 2-10 5 [45l.[60l|. 
For further evidence for the weak iia-dependence of the center 
temperature we also refer to table ITT1 



perature shift T c — T m and take full notice of violating 
OB. 3, which in contrast assumes constant dp/dT. We 
remind that this cannot be taken into account in the ex- 
tended BL theory presented in section [TTJ while DNS can 
well include it. 
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A. Numerical simulation approach 

To handle the numerical effort we restrict ourselves to 
incompressible and even two dimensional flow. The equa- 
tions governing non-Oberbeck-Boussinesq convection in 
incompressible fluid flow are the incompressibility condi- 
tion 

d lUl = , (14) 
the Navier-Stokes equation 

Pm(dtUi + UjdjUi) = -dip + dj(rj(djUi + diUj)) 

+gPm (1 - P/Pm) (15) 

and the heat-transfer equation 

{d t T + UjdjT) = dj(AdjT). (16) 

Here, <5^3 is the Kroneker symbol. The density is as- 
sumed to be constant and its value p m is fixed at that 



FIG. 12: Temperature dependences of the material properties 
v, k, p in the relevant T-range. dp/dT is also displayed, (a) 
Ethane around the temperature T m — 27° C adapted from 
[13. The pressure is P/P, = 1.062. (b) and (c): The fluid 
properties for water (b) and glycerol (c), respectively, around 
the temperature T m — 40° C (which has been studied in [45| 
and Note the significant variation of the density p with 

T in the case of ethane, as compared to the two liquids. In 
ethane, p/ ' p m — 1 various from about 0.07 for T — T m — —5K 
to about —0.1 for T — T m — +5K, whereas in glycerol prho m 
is practically independent of T. The strong variation of p for 
ethane follows from the large dp/dT (also shown). 



of the temperature T m , except in the buoyancy term, 
where the full nonlinear temperature dependence of p(T) 
is implemented. The dynamic viscosity r](T) and the 
heat conductivity A(T) are also both temperature and 
thus space dependent. The isobaric specific heat capac- 
ity cp is assumed to be constant, its value being c p . m 



11 



(in contrast to real ethane). The experimentally known 
temperature dependences of p, 77, A and the values of the 
parameters p m , c p ^ m for ethane are given in J4?J and, for 
better reference, are reported in Table U in the specific 
form implemented in our DNS. 

For consistency with the experimental measurements 
and with the BL theoretical analyses presented above for 
liquid-like ethane, we chose the arithmetic mean temper- 
ature to be T m = 27°C and the pressure as P/P* = 1.062. 



B. Numerical results: T c shift in liquid-like ethane 

From figure im we can conclude that the DNS captures 
the experimental measurements of the center tempera- 
ture shift T c — T rn as a function of A quite reasonably. 
The quality of the agreement with the available ethane 
data is similar to the one we have observed for glycerol 
[45l | and for water [6(|. This also serves as a further val- 
idation of our numerical approach. 

We note that for water and for glycerol the T c shift ob- 
tained by the extended BL theory [43 ] is nearly the same 
as calculated by DNS (see [45|,[6(|). In contrast, for the 
liquid-like ethane, the extended and even compressible 
BL theory only provides the right trend in the shift, but 
can not capture its amplitude (see again Fig. [TT|) . This 
observation supports our guess on the relevance of the 
nonlinear T-dependence of p(T) and thus of buoyancy. 
This latter is fully included in the DNS, as described in 
the previous subsection, while in BL theory it cannot be 
taken care of. 

One of the advantages of the DNS as compared to real 
experiments is that the material properties are easily and 
independently tunable. Therefore, the dynamics of hypo- 
thetical ethane-like liquids can also be addressed. In the 
next subsection we shall see how this approach is useful 
in understanding the effects of the temperature depen- 
dence of the various material properties on the center 
temperature shift. 
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C. T c shift in hypothetical fluids 

To obtain more insight into the physical origin of the 
non-Oberbeck-Boussinesq T c -shift, we consider NOB cor- 
rections for hypothetical ethane-like fluids in which at 
least one of the temperature dependences of k(T), v(T), 
$(T) is switched off, and fixed at the OB values K m , 
c m , (3 m . The quantity $(T), defined in eq. 0, is useful 
for the classification of the hypothetical fluids discussed 
in the following sections, but (3(T) is not explicitly in- 
troduced into the DNS, in which the density difference 
p(T) — p m is taken instead, see Eq. (TTS]) and Tab. U Fi- 
nally, we remind that /3(T) = j3 m , the usual thermal 
expansion coefficient, if the fluid density is a linear func- 
tion of the temperature around T m (i.e., if the conditions 
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FIG. 13: The normalized center temperature shift 
(T c — T m )/A versus the temperature difference A for 
several hypothetical liquids. We consider the following six 
hypothetical liquids: 0(T), n(T),u(T)), (/3 m , n(T),u(T)), 
0(T),K(T),u m ), (p m ,n(T),v m ), 0(T), Km ,u(T)) and 
(/3 m , n m ,v(T)). In particular in each panel we compare two 
cases, which differ only by their buoyancy's T-dependence, 
i.e., (3(T) instead of /3 m . The symbols indicate the simu- 
lation results at Rayleigh number Ra = 10 6 . The circles 
(O) correspond to the cases in which the full temperature 
dependence of the buoyancy (T — T m ) j3(T) is taken into 
account, while the triangles (A) represent the cases where 
only the linear temperature dependence (T — T m )/3 m is 
considered. The solid line shows the prediction of boundary 
layer theory with an incompressible flow assumption [43| ]. 
The dashed line stems from the solution with no convective 
flow d z (n(T)d z T) = with boundary conditions T = TL at 
z = and T = T t at z = L. See also Table HU 
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TABLE I: Expansion coefficients of material properties of ethane around the temperature T m = 27°C adapted from [47| . The 
pressure normalized by its critical value is P/P, = 1.062. The effective kinematic viscosity, the effective thermal diffusivity, and 
the buoyancy are written in a polynomial form as v(T) = rj(T)/p m — ~^2 n - A n {T — T m ) n [m 2 /s], k(T) = A(T)/(p m c p , m ) = 
Z) n =o Bn ( T ~ Tm ) n [ m2 / s l. and ff(l - p(T) / p m )=Y, n =i Cn ( T ~ Tm )" [ m /s 2 T, respectively. Using the leading coefficient for the 
buoyancy force, we can write the Rayleigh number as Ra = i3mL 3 A/(i/ m Km), where j3 m — Ci, v m = Ao and K m = Bo, which 
coincides with the usual OB definition. The polynomial expensions for (3{T) and $(T) are g/3(T) = X^n=i C n {T — T m ) n_1 and 
gP(T) = ( Pm /p(T)) ■ £ n=1 nC n (T - T m ) n -\ 



OB. 3 holds). In that case the thermal expansion coeffi- 
cient P{T) = /3 m /[l - (T - T m )/3 m ] still depends on T 
unless /?A < 1. 

For convenience hereafter we will call the two classes of 
artificial fluids, based respectively on the full non-linear 
NOB buoyancy force and on the linear OB approximation 
as defined by OB. 3, as the /3(T)-fluids and /3 m -fluids. 

In figure [13] we present the DNS results of the normal- 
ized temperature shift (T c — T m )/A for several types of 
hypothetical fluids. The numeral values are given in ta- 
bleflTJ One clearly observes in the figures and in the table 
that the temperature dependence of the thermal expan- 
sion function (5(T) apparently is relevant for the shift of 
T c . What can also be noticed in particular from table ITT1 
is that the corrections of the center temperature originat- 
ing from the temperature dependence of either v, k, or 
f3 are approximately additive (i.e., add "linearly"): E.g., 
the center temperature corrections of the fluids with (/3 m , 
K m , v{T)) and with (/3 m , k(T), v m ) add to that of the 
fluid with (f3 m , k(T), v(T)), etc. Note that this additivity 
is in contrast to what had been found within the extended 
BL theory of ref. [45| where the full compressibility of 
the density had been taken into consideration. Instead, 
in DNS we have restricted to incompressible flow. 

1. P{T)- and (3 m - fluids and extended BL 

In figure [13] we also compare fluids of (3{T)- and fi m - 
type with equal properties of thermal conductivity and 
kinematic viscosity on the different panels. For conve- 
nient comparison the predictions by the extended BL the- 
ory under the assumption of fluid incompressibility [43| 
are also shown for the hypothetical fluids. Note that in 
such case the chosen type of buoyancy force needs not be 
specified, because BL theory does not capture it: For the 
BL theory (3{T)- and /3 m -ethane are not distinguishable 
with respect to buoyancy. As a reference, the T c -shift 
value in the purely conductive case [ui = everywhere) 



is also reported. The T c shift in this case is not zero for 
fluids with k — k(T) as a result of the solution of the 
heat conduction equation k d%T + ^(d z T) 2 = with 

We now discuss our main findings. The change in T c 
for the hypothetical liquids with (3 m relative to the ones 
with P(T) is comparable for given A. More importantly, 
we find that the deviation (T c — T m )/A calculated with 
DNS is well captured by the BL theory for the cases of 
artificial /? m -fluids, i.e., for fluids where p(T) is assumed 
to be a linear function of T. In contrast, BL theory is 
always far from the T c shift obtained for the more real 
/3(T)-ethane. This indicates that the extended BL theory 
well captures the NOB effect once the assumption (OB. 3) 
is satisfied, even if (OB.l) and (OB. 2) are violated, but 
it does not correctly describe the NOB effects, if (OB. 3) 
is violated. Furthermore, the DNS results reveal that the 
T c -shift is always enhanced if a temperature dependence 
of the thermal expansion function (3(T) determines the 
buoyancy, i.e., if p(T) depends nonlinearly on T. 



2. Mirror transformation ${T) -» 0(2T m - T) 

To quantitatively appreciate the effect of the temper- 
ature dependence in each material property individually 
also at different Rayleigh numbers, Ra = 10 6 -10 8 , we list 
the T c shifts for several hypothetical ethane-like fluids in 
Tab. [IT] Since here our primary concern is the influence 
of the thermal expansion function, besides the (3(T)- and 
/3 m -fluids a new class of hypothetical fluids is introduced. 
We consider (3(2T m — T)-fluids obtained by the mirror 
transformation /3(T) — > (3(2T m — T). As schematically 
shown in figure 1141 this transformation reverts the non- 
linearity in the buoyancy force with respect to T — T m . 
The comparison between the cases (J3(T), re m , u m ) and 
((3(2T m - T) ) - cases 4 and 8 in table [IT] - shows 

that the effect of the mirror transformation on the out- 
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case 


$ 




K 


V 


100(T C - T m )/A 
at Ra = 10 6 


100 {T c — T m )/A 
at Ra = 10 7 


100(T C - T m )/A 
at Ra = 10 8 


1 (NOB) /3(T) 




k(T) 


u(T) 


1.3003 ± 0.0369 


1.3879 ± 0.0527 


1.0788 ± 0.0362 


2 


$(T) 




k(T) 




0.2699 ± 0.0356 


0.2174 ±0.0515 


0.1251 ±0.0414 


3 


Kt) 




ftm 


u{T) 


2.4283 ± 0.0361 


2.4351 ± 0.0496 


2.4370 ± 0.0462 


4 


P{T) 




ft m 


Vm 


1.4805 ± 0.0355 


1.4320 ± 0.0735 


1.3796 ± 0.0534 


5 


P(2T m 


-T) 


k{T) 


v{T) 


-1.5953 ±0.0361 


-1.7946 ±0.0506 


-1.5868 ±0.0320 


6 


P(2T m 


-T) 


k(T) 


v m 


-2.6382 ± 0.0357 


-2.7458 ± 0.0496 


-2.5516 ± 0.0496 


7 


P(2T m 


-T) 


ft m 


u{T) 


-0.4738 ± 0.0375 


-0.4578 ± 0.0580 


-0.3358 ± 0.0354 


8 


P(2T m 


-T) 


ft rn 


v m 


-1.4878 ±0.0366 


-1.5358 ±0.0504 


-1.4013 ±0.0363 


9 


(3 m 




K(T) 


v{T) 


-0.1983 ±0.0361 


-0.2043 ± 0.0467 


-0.2691 ± 0.0587 


10 


Pm 




«(T) 


Vm 


-1.2369 ±0.0363 


-1.1883 ±0.0611 


-1.1851 ±0.0369 


11 


Pm 




ft-rri 


u{T) 


0.9852 ± 0.0364 


1.0834 ± 0.0576 


0.9916 ± 0.0259 


12 (OB) 


(3-m 




ft-rn 




0.0171 ±0.0381 


-0.0484 ± 0.0548 


0.0271 ± 0.0387 



TABLE II: DNS results for the center temperature shift T c — T m normalized by the temperature difference A = 10K 
for several hypothetical fluids. The effective thermal expansion function is given by f3(T) = g" 1 ^2 —i Cn(T — T m )" _1 
[1/K]. Using the expansion coefficients C„ listed in Table U we write the buoyancy for the case of /3(2T — T m ) as 
<7(l-p(T)/p m )=£ n=1 (-l) n+1 a i (T-T m r [m/s 2 ], and in the m case as g(l - p{T)'/ p m )=C\(T - T m ) [m/s 2 ]. Although the 
center temperature shift for the OB case (case 12) should be essentially zero due to the top-bottom symmetry, the mean value 
determined from the DNS result is non-zero since the sampling time for taking the statistics is finite. Note that measurements 
of the temperature shift in the OB case are all, within statistical uncertainty, compatible with zero. 





flCO / 






j3(2T m -T)\ 





-A : ^ T 

T T T 



FIG. 14: Schematic plot of the mirror transformation f3(T) 
P(2T m — T) of the thermal expansion coefficient /3. 



put parameter T c is to change the sign of (T c — T m )/A 
while preserving its modulus. Furthermore, the deviation 
of (T c - T m )/A for all /3 m -fluids (cases 9,10,11) relative 
to the /3(T)-fluids (cases 1,2,3) is always positive, while 
it is always negative relative to the (3{2T rn — T)-fluids 
(cases 5,6,7). These features hold at all the studied Ra 
numbers. Therefore, we conclude that the shift T c — T m 
is sensitive to the sign of the slope of (3(T) or, equiva- 
lently, to the sign of the nonlinear term in the buoyancy 
factor 1 — p(T)/p m . More precisely speaking, the mirror 
transformation changes the signs of the even order co- 
efficients C*2 and C4 defined in Tab. |TJ Obviously Ci is 
the larger and thus the most relevant coefficient. These 
features are absent in water and in glycerol, because for 
those the temperature dependence of p(T) is much less 
pronounced. 



3. Test on the linearity of the T c — T m shift 

Looking at Tab. |TT] in more detail we find that all 
changes in (T c — T m )/ A, which stem from the nonlinear 
T-dependence of the buoyancy force, i.e., from the dif- 
ferences between (3 m and /3(T) (or (i{2T m — T)) but hav- 
ing the same temperature dependences of k and are 
comparable. To emphasize this feature, we look at the 
differences T C \^ T) - T c \ 0m (or -{T c \^ 2Tm _ T) - T c \p n J) 
and plot them versus A as well as versus Ra, see fig- 
ure [T3J We find a good collapse of the data onto a 
single curve for various temperature dependences of k 
and v. In particular, the comparison between the mid- 
dle and bottom panels of Fig. [T5] leads to the relation 
T "\$(T) ~ T c\p m = _ ( T 'c| / §( 2Tm _ T) - T c \p m ), which indi- 
cates that the T c change is dominated by the quadratic 
term in 1 — p(T)/p m , but is almost independent of k 
and v. This observation may be important for further 
attempts to improve the extended NOB BL theory. 

VI. CONCLUSIONS 

In this paper we have first presented in full detail the 
extension of boundary-layer theory to the case of com- 
pressible NOB fluids in a Rayleigh-Benard system. The 
theory predicts a deviation of the center temperature T c 
from the arithmetic mean temperature T m between the 
top and bottom plates, i.e., T c — T m ^ 0. 

Second, the theory has been tested against new exper- 
imental data for ethane near the critical point in its more 
liquid-like phase. Data come from a series of experiments 
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FIG. 15: Effect of the temperature dependence of the ther- 
mal expansion coefficient (3 on the shift of the center tem- 
perature for several hypothetical fluids. DNS results for the 
normalized temperature difference (T C \^( T \ — T c \p m )/A (or 

— (T c \p(2T -t) ~ T c \p m )/A) with the same temperature de- 
pendence of the thermal diffusivity k and the kinematic vis- 
cosity v are plotted; here T c \p, T -. denotes the center temper- 
ature with the full temperature dependence of the buoyancy 
g(l — p/pm) as given in Tabled and T c \/3 m denotes that with 
the linear temperature dependence g(l — p(T) / p m )=Ci(T — 
T m ) only. The top panel shows {T C \^ T ^ ~ T c \p m )/A versus 
A at fixed Rayleigh number Ra = 10 6 . The solid line shows 
the linear fit 1.473 • 10~ 3 K~ 1 x A for the case $(T),Km,v m 
(see Table|Il]|. The middle panel shows {T C \^ (T) - T c \p m )/A 
versus Ra at fixed temperature difference A = 10K. The 
bottom panel has the same parameters as the middle one, ex- 
cept showing — (T c \^ 2 t -t) ~ \p™ ) /A for the case of mirror 
transformed f3. 



in cylindrical cells of aspect ratio T = 0.5 and 1, reaching 
Ra numbers O(10 10 ). The experimental measurements 
at T m = 27°C, P/P* = 1.062 (and T = 1) have been cho- 
sen for comparison with those at T m — 40° C. Contrary 
to the good agreement observed for the case of gas-like 
ethane [431 ] , the BL theory here gives much smaller values 
of the center temperature shift as experiment. 

Third, direct numerical simulations DNS, based on T- 
dependent material parameters but still within the in- 
compressible approximation and a two-dimensional do- 
main have been performed to get more insight into the 
observed discrepancy between experiment and extended 
BL theory. The DNS results provide a satisfactory agree- 
ment with experiment both in the gas-like as well as in 
the liquid-like cases. Several hypothetical ethane-like flu- 
ids have been investigated too. Our analysis shows that 
the extended BL predictions fail whenever the non-linear 
temperature dependence of the density p(T) is imple- 
mented in the numerical simulations. Furthermore, if the 
dependence of f3(T) on T dominates the NOB effects, the 
sign of the linear term in the effective expansion function 
j3{T) is responsible for the sign of the variations of T c 
as compared to the OB value T c = T m . But even if the 
difference (3{T) — [3 m is a non linear function of the tem- 
perature, the tiny T c shift detected in our simulations is 
proportional to T c \^, T s - T c \p m . 

Our results can be summarized in the following physi- 
cal picture of the surprising phenomenon that the center 
temperature T c in liquid- like ethane near its critical point 
increases, the center becomes warmer than the arithmeric 
mean temperature T m , while in gas-like ethane the center 
gets colder, T c is less than T m . Namely, in the liquid-like 
case the buoyancy cx (T—T m ) j3(T) is larger at the bottom 
and smaller at the top, supporting the uprising warmer 
plumes more than the down coming colder top plumes. 
This brings predominantly hotter material into the bulk. 
For gas-like ethane the buoyancy is larger at the cooler 
top, which accelerates the downgoing cold plumes with 
more preference than the uprising warm plumes from 
the bottom, that experience a weaker buoyancy. This 
in turn brings more cooler material into the bulk, lead- 
ing to T c < T m . It is the sign of the slope of $ (or of (3), 
which is the relevant quantity. 

The more general lesson which can be drawn from this 
paper is that there is a plethora of origins of NOB correc- 
tions. Which one dominates can only be said by having 
a detailed look at the temperature dependence of all ma- 
terial parameters. Both the extended BL theory and 2D 
DNS are useful tools to judge which temperature depen- 
dence is the most relevant one or whether they all matter, 
as we now have often seen. 
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APPENDIX A: BOUNDARY-LAYER EQUATIONS 

1. Viscous boundary-layer 

Consider two-dimensional flow over a flat plate, such 
that the main wind velocity U does not depend on x and 
linXj—voo u x (x, z) — U. Then, the s-momentum near the 
wall is governed by [50| 



p\u x 



du x 
dx 



du x 
dz 



d 2 u x d 2 u x 



dx 2 dz 2 
dr\ ( du x du z 
dz 1 dz dx 



(Al) 



d_ 

dx 



n 



+ ( 



du x 
dx 



du z 
dz 



where r\ is the dynamic shear viscosity and £ the volume 
expansion viscosity. 

To appraise the dominant structure of equation (|Aip at 
sufficiently large Reynolds numbers, we follow Prandtl's 
scaling: 



are of order 1; the remaining terms are of order 1/Re.* 

Therefore, at large Re, the dominant part of the x- 
momentum equation is given by 



du x 
dx 



du x 
u z — 
dz 



d 2 u 



dr) du x 
dz 2 dz dz 



1. Thermal boundary-layer 

Consider again two-dimensional, subsonic, and steady 
flow over a flat plate. In the framework of boundary- 
layer theory, energy conservation leads to the following 
equation for the entropy per mass s [4l[ : 



rp ds 8 f dT 

u x pT V u z pT — = ^-<A — 

dx dz dz dz 



(A7) 



Letting s = s(T, p), the entropy gradient (say, its d z s 
component) can be expressed as 



ds 

dr: 



ds\ dT , fds\ dp 
dp) T dz 



dTj p dz 



The first contribution is directly associated with the iso- 
choric specific heat (per mass) of the gas 



dT 



c P 

c v = — ■ 

7 



L x, 
L 



/Re 
U u x , 
U 



z, 



Re 



(A2) 
(A3) 
(A4) 
(A5) 



The second contribution follows from a Maxwell relation, 



ds 

dp j T 



1 fdP 



P 2 \dT y p 



1 cp 7 — 1 



Thus, the left-hand side of equation (|A7|) can be written 
as 



where L denotes the typical length of the plate in flow di- 
rection and Re = LU p m /Vm the Reynolds number. The 
index m indicates that the fluid properties are evaluated 
at a thermodynamic reference state (T m , P m , p m ), which 
is adopted in the nondimensionalization of p = p m p and 
t\ = r) m fj. Then, substituting into (fATj) . one 

finds 



» du 
p< u x 



dx 



du x 
dz 



Re 9a; 2 
du x 
dz 



dr\ 
dz 
J_d_ 
Re dx 



d~z 2 
1 du. 
Re dx 





du x 


du z 




1% 


dx 


dz 


} 



ds ds pep / dT dT 
Ux pT h u z pT — — = <u x — +u z — 

dx dz 7 dx dz 

c P 7 - 1 f dp dp 
— < u m — + u 

7 



(3 \ x dx z dz 



Finally, using the continuity equation u x d x p + u z d z p 
-p{d x u x + d z u z }, one finds 



dT dT 7 - 1 f du x du z \ 7 d ( , dT 



dx dz \ dx dz J pep dz \ dz 

Note that the limiting case of liquids (namely 7 = 1) is 
fully accounted by this equation. 



where rj = f) ^ + ^ • Clearly, all terms on the left-hand 

side of equation (1A6|) are of order unity. However, this 
is not the case on the right-hand side of (|A6|) : only the 
viscous contributions involving transverse gradients of u x 



t Note that the term involving V ■ u : 



+ is of order I/Re, 
as long as £ and r) are of the same order of magnitude. Indeed, 
(~/t] is of order unity if acoustical effects (Ma 1) and chemical 
reactions do not take place in the fluid. For situations in which 
C > T), see [U Section 81]. 
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APPENDIX B: SELF-SIMILARITY ANSATZ 
1. Viscous boundary-layer 

In the stream- function representation (j8]) , the longitu- 
dinal velocity u x is expressed as 



insertion of (TB2]) and (|B3j) into (jB6|) leads to 
pi 2 d ( du x 



■q m U c dz \ dz 



77—^ \ = fj^f' 



Id ~ 

p 



p 9z 

and its spatial derivatives are 









P . 






r 






p 






p 



"'" |{P*"-P'*'}, (Bl) 



du x 

~~dx~ ~ 2£ 2 P 2 

<9z 2 



Therefore, by substituting (|B5|) and (|B8j) into ([SJ) one 
finally obtains equation (jlip . 

2. Thermal boundary- layer 



4 p 2 

e p 3 1 



p 2 *"'-2pp'*" 



In the same spirit as above, inserting (|B 1[) . (|B4[) , and 
A™ A 



— (a — 

dz I 



G 



{ao" + a'o'} 



+ [2(p') 2 - P~p"\ • (B3) i nto equation © one finds 



Likewise, from equation ©, the transverse velocity u z 
reads 

pox 24 \ p p \ 



e" + A'e' 



with 



du z v m 1 r - - w - , ~ , 



(B4) 



Thus, the advective contributions in Prandtl's equation 
([5]) can be written as 



du x du x 
P i u x — h Wz-^— 

02; az 



1 ?7m 



2 pP c 



— * * ' \ . (B5) 



Now consider the viscous contributions: 

d ( du x \ d 2 u x d-q du x 

v^rr = V^OT + — -^r- (B6) 



dz 1 dz 



dz 2 dz dz 



Since 



dri t) m J / di) \ ~, ( df}\ 



(B7) 



2^ p 2P C (3 p 2 pc P 
where /3 = /3 A and Pr = (z/ m p m cp yTn )/ A m . Thus 

A0'' + J±^p r * + A'le'-^^Pr* = 0. 
I27 J 2 7 p 

Here, substituting p' by (fT2")l one finally obtains 



A9" + \ -c P Pr * + A' } Q' 



Note that the limiting case of liquids (namely 7 = 1) is 
fully accounted by this equation. 



APPENDIX C: NUMERICAL RESULTS ON 
NUSSELT NUMBERS IN REAL AND 
HYPOTHETICAL ETHANE FLUIDS 

For completeness, in table UlTl the NOB corrections in 
the Nusselt number are given, resulting from the numeri- 
cal simulations of real and hypothetical ethane. The cor- 
responding NOB corrections of the center temperature 
had already been shown in table [TTJ 
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TABLE III: Simulation results of the relative deviation of the Nusselt number Nu/Nuob — 1 at the temperature difference 
of A = 10K for several hypothetical fluids. At the Rayleigh numbers of Ra = 10 6 , 10 7 , and 10 s the OB Nusselt numbers 
correspond to Nuob = 6.53, 12.43, and 25.12, respectively. The notations for the material properties are the same as in Table 
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